


x 



On an integrable discretization 
of the modified Korteweg-de Vries equation 



Yuri B. SUMS 



Centre for Complex Systems and Visualization, University of Bremen, 
Universitatsallee 29, 28359 Bremen, Germany 
e-mail: suris @ cevis.uni-bremen.de 



> 

m 

o 
o 

(N 

o 
l> 

a^ , 

Abstract. We find time discretizations for the two "second flows" of the Ablowitz-Ladik 
hierachy. These discretizations are described by local equations of motion, as opposed to the 
previously known ones, due to Taha and Ablowitz. Certain superpositions of our maps allow 
a one-field reduction and serve therefore as valid space-time discretizations of the modified 
(y:j . Korteweg-de Vries equation. We expect the performance of these discretizations to be much 

better then that of the Taha-Ablowitz scheme. The way of finding interpolating Hamiltonians 
for our maps is also indicated, as well as the solution of an initial value problem in terms of 



Lj ■ matrix factorizations. 



1 Introduction 

Already in the early period of the soliton theory it was realised that by the numerical simu- 
lation of soliton equations it is highly desirable that the difference schemes inherit the inte- 
grability property [0], P], [Q. These papers contained a full (space and time) discretization 
of such soliton equations as the nonlinear Schrodinger equation (NLS), Korteweg-de Vries 
equation (KdV) and the modified Korteweg-de Vries equation (niKdV). This was done 
in two steps. On the first step the space discretization was performed |T]]. To this end 
the authors discretized the auxiliary linear problem related to the corresponding soliton 
equation. Concerning as an example the mKdV equation, 

qt = qxxx =F Qqxq^ (1-1) 

the associated linear problem is the Zakharov-Shabat problem 

( K q 
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Its discretization chosen in |]l| is 

/A qk 
^fc+i = I ^fc, (1.3) 

V Mk A-i 

and the space discretization of (|1.1| ) found in 0] is 

qk = (1 =F ql) (qk+2 - 2gfc+i + 2qk-i - qk-2 

=F qk+i{qk+2qk+i + qk+iqk) ± qk-i{qkqk-i + qk-iqk-2)) (1-4) 



(to perform the corresponding continuous limit, one has to assume in ( |1.4| ) qk = eq{ke), 
to rescale the time t \—> t/(2e^), and then to send e — > 0). 

The second step of discretization - the time discretization - was performed in (for 
NLS) and in 0] (for KdV and mKdV). The approach to this second step of the discretiza- 
tion process was fundamentally different: the linear problem ( |1.3D was not modified any 
more, and only a suitable choice of a (discrete)-time evolution of the wave function \l'fc was 
imposed. The outcome of [Q was an excellent numerical scheme for the mKdV, though 
not free of some drawbacks, the main one being the non-locality (which means that qk, 
the discrete time update of qk, depends explicitly on all g^-'s and q/s with j < k). This 
property is unpleasant from the aesthetical point of view as well as from the practical one 
(since it implies a large amount of computations by numerical realization). 



In a more modern language, time discretizations in [^, were sought in the same 
hierarchies to which the underlying continuous time systems belong. Recently this ap- 
proach was pushed forward as a systematic procedure of obtaining integrable discretiza- 
tions @]-[^. The results of [§ were re-considered in |§], where a significant amendment 
was achieved. Namely, the non-locality of the schemes in p[ for NLS was overcome. In 



the present paper we do an analogous work for the time discretization of (|1.4|) from 



2 Ablowitz— Ladik hierarchy 

To deal with the equation ( p..l|) in a slightly more systematic way, one has to consider it 
as a particular case of the following, more general system: 

qt = Qxxx - Qqxqr, rt = r^^^ - Gqr^r (2.1) 

under the reduction 

r = ±q. (2.2) 

Analogously, the space discretization ( |1.4|) arises by the same reduction from the more 
general system 

Qk = (1 - Qkrk) {qk+2 - 2gfc+i + 2qk-i - qk-2 
-qk+i{qk+2rk+i + qk+irk + qkTk-i) 

+qk~i{qk-2rk-i + qk-iTk + qkTk+i)) 

(2.3) 
Tk = (1 - qkTk) (rfc+2 - 2rfc+i + 2rk-i - rk~2 
-rk+i{qk+irk+2 + qkTk+i + qk~irk) 

+rk-i{qk-irk~2 + qkTk-i + qk+iTk)) 

Below we consider this system either on an infinite lattice {k G Z) under the boundary 
conditions of a rapid decay (l^fcl, I'^fcl — > as A; ^ ±oo), or on a finite lattice {I < k < N) 
under the periodic boundary conditions (go = qN, f^o = rjy, qN+i = qi, ^n+i = ^i). In any 
case we denote by q (r) the (infinite- or finite-dimensional) vector with the components 
qk (resp. r/,). 

From the modern point of view, the system ( p.3| ) is a representative of a whole hierarchy 
of commuting Hamiltonian flows, the Ablowitz-Ladik hierarchy. An object of the principle 
importance in the description of this hierarchy is the 2x2 Lax matrix 

A qk 
Lfc = Lfc(q,r)=| | (2.4) 

Tk A-i 



depending on the variables q, r and on the additional (spectral) parameter A. Each flow 
of the hierarchy allows a commutation representation (semi-discrete version of a zero- 
curvature representation) 

Lk = Mk+iLk - LkMu (2.5) 

with some 2x2 matrix M^ = Mfc(q, r, A). 

The Hamiltonians of the commuting flows are the coefficients in the Laurent expansion 
of the trace tr T/v(q, r, A) where T/v is the monodromy matrix 



■N 



Ln ■ Ln-1 ■ 



. . . ■ L2 ■ Li, 



supplied by the function 



N 



ifo(q, r) = logdetTAT = ^ log(l - g^rfc). 

fc=i 

The Poisson bracket on the phase space is given by 

{lk, rj} = (1 - qkrk)Sjk, {qk, Qj} = {rk, rj} = 0. 



(2.6) 



(2.7) 



(2i 



The involutivity of all integrals of motion follows from the fundamental r-matrix 
relation: 

{L(A) f L(/x)} = [L(A) ® Lifi), p(A, /.)], (2.9) 

where 

\ 

T^ -h 



' 2 A2-|t2 



p(A,/i) 



v 



2 A2-Ai2 

X/j. _1 

A2-Ai2 2 

"^ iA!±m! 
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2A2-u2 / 



We refer the reader to |^ for the general formula expressing the matrix M^ in terms 
of the corresponding Hamiltonian function -ff(q, r) and the r-matrix p(A, /x). 

It is easy to see that the following functions belong to the involutive family generated 
by tr Tn: 

N N 

Hiici, r) = J2 Qk+iTk, ^-i(q, r) = Zl Qkrk+i, (2.11) 

fc=l k=l 



N 



N 



N 



^2(q,r) 
^-2(q,r) 



J2 1k+2rk - Y. Qk+2qk+irk+irk ~ 7;Y1 1k+i 



2 2 
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k=l 

N 



k=l 

N 



k=l 

N 



Y. 1krk+2 - Y Qkqk+irk+irk+2 - ^ H lM+v 

k=l k=l ^ k=\ 



f2.12l 
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The corresponding Hamiltonian flows are described by the differential equations 



qk = g/c+i(l -qkrk) 
Tk = -rfc-i(l - QkTk) 



-^-1 



qk = Qk-iC^ - Qkrk) 
Tk = -rk+i{l -QkTk) 



To 



T. 



-(^fc-2 - rk-i{qk-irk-2 + QkTk-i + Qk+iTk 
-- (qk~2 - qk-i{qk-2rk-i + Qk-iTk + qkrk+ 



1 - Qkrk) 
(1 - qkVk] 

1 - Qkrk) 



(2.14) 
(2.15) 
(2.16) 

(2.17) 



Qk 

Tk = -{rk+2 - rk+i{qk+irk+2 + qkTk+i + qk-iTk)) (1 - qkTk 

The flow ( ^.31 ) is an obvious superposition of these more fundamental and simple flows, 
namely 

T2{t) O ^i(-2t) O ^_2(-t) O ^_i(2t) 

In what follows we will concentrate on the flows 

-^2(t)oJ-i(-ct) (2.18) 

^_2(t)o^_i(-ct), (2.19) 

having in mind that in the application to the mKdV case the value c = 2 is of interest. 
The matrix Mk corresponding to the flow ( p.l8| ) is given by the formula 

/ A^ - X'At^ - Af \\k + Ai3« \ 

\ x^Tk-i + xc'i^ x^kTk-i ) 



r(2) 



M):\c) 



(2.20) 



with 



^(2) 



k 

(1) 
k 

(1) 



c + qkTk-x 

qk+i - qk{c + qk+iTk + qkTk-i) 

rk-2 - rk-iic + qk-irk-2 + Qk^k-i) 

qk+iTk-i + qkrk-2 - qkrk-i{c + qk+iTk + qkTk-i + qk-irk~2) 

Similarly, the matrix Mk corresponding to the flow ( p.l9|) is given by 

/ -X'^qk^iTk -XSk-i-X-^B[7^^ \ 



C 



AT 



Mt\c) 



V -A -"rk 



x-'ct' 



-x-^ + x-^vi-'^ + vT 



(2.21) 



with 



V 
B 



(-2) 
k 

(-1) 



c 



(-1) 



V 



(0) 



c + qk-iTk 

qk-2 - Qk^iic + qk-2rk-i + Qk-iTk) 

Tk+i - rk{c + qkTk+i + Qk^iTk) 

qk-2rk + Qk^iTk+i - Qk-irkic + gfc-2^fc-i + Qk-irk + Qk^k+i)- 



Obviously, one has: 

Mf ^ (c) = Mf ) (0) - cM^^^ , Mi-'^ (c) = m(-^) (0) - cUi-'^ , 
where the matrices 



M, 



M, 



(-1) 



(1) 



/ 



A^ - QkTk-l Mk 

Xvk-i 

-\-\k- 



\ -A Vfc -A 2 + (?fc-irfe 



correspond to the flows J-'±i. 



The matrix Mk for the system (|2.3| ) is equal to 



(2.22) 



(2.23) 



3 General remarks about the time discretization 



In 1^ Taha and Ablowitz constructed a time discretizations of the system ( |2.3| ), thus 
achieving a full discretization of the system (|2.1| ). The basic feature of the time discretiza- 
tion in P] is following: it admits a discrete analog of the zero-curvature representation, 



LkVk — Vk+iLk 



(3.1) 



with the same matrix Lk as the underlying continuous time system. (In ( p.l|) and below 
we use the tilde to denote the /i-shift in the discrete time /iZ). In a more modern language, 
the maps generated by the discretizations in ||^ belong to the same integrable hierarchy 
as the continuous time system (|2.3|). 

We shall not need concrete expressions for the entries of the matrix Vk and the cor- 
responding evolution equations obtained in [Q. However, the details of the derivation 
seem to be never published, therefore we present them (in a slightly amended form) in 



the Appendix. Here we restrict ourselves with some general remarks following from the 
work of Taha and Ablowitz and necessary for the following presentation. 

They considered difference equations allowing a commutation representation ( |3.1j ) with 
the matrix L^ given by ( p.4|) and assumed that the entries of the matrix 




Vu=\ I (3.2) 

have the following A-dependence: 



Au{X) - 


- 1 + h\'A^^ + /.AMf + hAf + h\~'A[~'^ + hX-'A[''^ 


Du{\) - 


- 1 + hX'D^^ + hX^Df^ + hD^^ + hX-^D^^^^ + hX'^D^- 


BkiX) -- 


-- hX'B'^^ + hXB'^^ + hX-'Bt'^ + hX~'Bi-'^ 


Cu{\) - 


= /lA^cf ) + hXCl^^ + hX-'ct'^ + hX--'ct'\ 



They showed that each such difference equation may be completely characterized (in 
the case of rapidly decaying boundary conditions) by the limit values 



a(^) = lim A'j^\ 6^^^= lim D^^^ (j = 4, 2, 0, -2, -4) (3.3) 

fc— >itoo fc^itoo 

Under the condition 



a 



(i) 



5^-^^ (j= 4, 2, 0,-2, -4) (3.4) 



the corresponding difference equation allows the reduction 

r = ±q. (3.5) 

The last statement is easy to see directly. Indeed, the discrete zero-curvature equation 
( p.l| ) is equivalent to the following four equations: 

^fc+i(A)-Afc(A) = X-\quCk{X)-nBk+i{X)) (3.6) 

Dfc+i(A)-Dfc(A) = A(ffc5fe(A)-gfeCfc+i(A)) (3.7) 

A-i5fc+i(A)-Afifc(A) = gfcA(A)-gfeA,+i(A) (3.8) 

XCk+i{X)-X-^Ck{X) = 7kAk{X)-rkDk+i{X) (3.9) 

Obviously, these equations allow the reduction 

r = ±q, A{X) = D{X-^), B{X) = ±C{X-^). 



But in this reduction, obviously, the condition ( p.4|) is satisfied. Since the difference 
equation is completely characterised by the quantities (^.31), the condition ( ^.41 ) is also 
sufficient for the above reduction to be admissible. 

Equating coefficients by different powers of A in (|3.6| )- (|3.9| ), one obtains 20 equations. 
In 1^, Taha and Ablowitz assumed a special role of evolution equations for g^, fk to 
two of these equations, using the other 18 to determine the 18 unknown functions A^\ 
ni^^ {j = 4, 2, 0, -2, -4), and b''^\ ci^^ (j = 3, 1, -1, -3). This way of dealing with the 
problem (inherited from the continuous time case) resulted in highly nonlocal expressions. 
This feature of the resulting difference scheme made its numerical realization extremely 
time consuming. Nevertheless, even despite the drawback of non-locality, this difference 
scheme proved to be the best among the numerical methods tested in . The reason for 
this lies undoubtedly in the integrable nature of this scheme. 

The goal of the present paper is to demonstrate how the feature of nonlocality may 
be overcome. This is achieved with the help of two basic ideas. First, we factorize the 
difference scheme into the product of several simpler ones, corresponding to the funda- 
mental and simple ffows of the Ablowitz-Ladik hierarchy. Second, by dealing with these 
simpler schemes we consider all the 20 equations on an equal footing. This allows us to 
derive simple local schemes, which we expect to exceed the original one due to Taha and 
Ablowitz. 



4 Local discretizations of the flows J^±2 

In the present Section we introduce two maps which serve as time discretizations of the 
ffows i ^JBi ), ( ^191) . 



Theorem 1. Consider the map 7^(/;,;c) : (q, r) \-^ (q? r) defined by the following 
equations of motion: 



T2{h-c) : < 



(?fc - Qk)/h 



{rk - rk)/h 



{qk+2 - Qk+iPk+i) (1 - qkrk) 

(1 - hqu+iTk-i) 

(4.1) 
{rk-2 - Tk-lPk) (1 - qkTk) 



(1 - hqk+irk-i) 

where the local quantity P^ = P{qk+i, qk, Qk-i, ^k, rk-i, ^fc-2; h; c) is defined by 

[1 - hqkqk-irk-irk-2) (1 - hqk+iqkhrk-i) 



{Pk - hqk.A-2] ^^ _ ^^^^^^^^ ^^ _ ^^^^^^^_^^ 

- [qk+iTk + qkTk-i + qk-irk-2 - hqk+iqkqk-irkrk-irk-2) = c, (4.2) 



so that 



Pk = c + Qk+irk + qkTk-i + qk-irk-2 + 0{h). 



These equations approximate the flow ( p.l8| ) and have the commutation representation 
^AI) with the matrix 



v; 



(2) 



^ l + h\^~ hX'^Af - hAf hX^qk + hXBl}^ \ 
^ hX^h-i + hXCi^^ l + hX^kh-i ) 



(4.3) 



Here 



.(2) 
^k 

4(0) 



Pk — qk+l^k — Q'fc-l^fc-2 

qk+i - qk{Pk - qk-irk-2) 
rk~2 - h-iiPk - qk+iTk) 
qk+ih-i + qkrk-2 - qkh-iPk 



(4.4) 
(4.5) 
(4.6) 
(4.7) 



Theorem 2. Consider the map T^2{h;c) : (q, r) t— > (q, r) defined by the following 
equations of motion: 



T^2{h]c): < 



{qk -qk)/h 



{rk - rk)/h 



{qk-2 - qk-iSk) (1 - qkTk) 

(1 + hqk-irk+i) 
{rk+2 - Tk+iSk+i) (1 - qkTk) 



{U 



(1 + hqk-irk+i) 

where the local quantity Sk = S{qk-2, qk-i, qk, ^k-i, rk, rk+i, h; c) is defined by 

(n ,f^~ \ (1 + hqk-2qk-irk-irk) (1 + hqk-iqkrkrk+i) 

^ ^ ^=-2 k+ij ^^ ^ hqk_2rk) (1 + hqu-irk+i) 

- {qkTk+i + qk-iTk + qk-2rk-i + hqk-2qk-iqkrk-irkrk+i) = c, (4.9) 

so that 

Sk = c + qkTk+i + qk-iTk + qk-2rk-i + 0{h). 

These equations approximate the flow ( p.l9|) and have the commutation representation 



( p.l[ ) with the matrix 

' 1 - hX-^k-iTk -hX~'^qk-\ - hX-^B^. 

^ -hX-^rk - hX-^ct^^ 1 - hX-^ + hX-^D['^^ + hof^ ) 



V, 



(-2) 



\ 



(4.10) 
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Here 

d[~^^ = Sk - qkTk+i - qk-2rk-i (4.11) 

bIt^^ = qk-2 - qk-i{Sk - qtrk+i) (4.12) 

Ck = Tk+i - rk{Sk - qk-2rk-i) (4.13) 

Df^ = qk-iTk+i + qk-2rk - qk-iTkSk (4.14) 

Proof. Since the both Theorems are proved similarly, we give here only the proof of 

the Theorem 2. Substituting the ansatz ( 4.101 ) for the matrix V^ into ( p.l|) , one sees 



that the following 6 equations have to be satisfied: 

^i+?-^i"^^ = qkrk+i - qk-ih (4.15) 

nn^-D^k' = ^.Cil'^-hBt^ (4.16) 

4;!^ = qk-i-qk{qkrk+i + Di-'^) (4.17) 

Ct'^ = rfc+i - rfc (g-fc.iFfe + 4;?) (4.18) 

qk-qk = hB[-^^ - hqkDf'^ (4.19) 

h-Tk = -hCi-\^ + hrkD'f'l, (4.20) 

As indicated above, the approach by Taha and Ablowitz to these equations would be as 
follows: consider the last 2 equations as the equations of motion, where the (non-local) 
expressions for the quantities Bf, \ Cl , -D^ follow directly from the first 4 equations. 
The crucial feature of our approach to the solution of these equations is that we use all 6 
equations to derive local expressions. In other words, we do not assume that the last two 
of them play a special role. 

Note first of all that due to ( |4.15| ), the equations ([4.17|) , ([4.18|) may be equivalently 



re-written as 

4;!^ = gfc_i-gfc(gfc_iFfc + Di;?) (4.21) 

Ct''^ = rk+1 - r-k {qkTk+i + Dir^^) (4.22) 

Now we introduce the auxiliary quantity Sk by the formula 

Sk = Dl ' + qkTk+i + qk-2rk-i- (4.23) 



which immediately gives (|4.11 ). With the help of this quantity we re- write the equations 



(g), (|0^ ) as (|;g) and QOBP , respectively. Substituting ( ^TT^ , ( [O^ ) into ( ^TT^ , we 
obtain the equation which may be written as 



i(o) n(o) 



Dll^ - Dl' = qkrk+2 - qk-2rk - qkTk+iSk+i = qk-iTkS, 



k 



-{qk - qk){rk+2 - rk+iSk+i) - (vk - rk){qk-2 - qk-iSk)- (4.24) 
9 



Similarly, substituting ([4.12|) , ( [4.13|) into (|4.19| ), ([4.20|) we obtain two equations which 
may be put down as 



^1 + hqk-^irk+i) {qk - Qk) = ^(gfc-2 - Qk-iSk) - hqk (of^ - qk-iVk+i 



(4.25) 

(1 + hqk-iTk+i) (ffe - Tk) = -h{rk+2 - Vk+iSk+i) + hvk (^i+i - qk-iTk+i) (4.26) 

Multiplying (f4.24|) by (1 + /igfc-i^fc+i) and using on the right-hand side the expressions 
( [4.25| ), (|4.26| ), we obtain an equation, which after some identical re-arrangements may be 
presented as 

^1 + /i/^i+i) (l + hqk-iTk+i + hqk-2rk - hqk-iVkSk] 



So, the quantity 



does not depend on k, and, taking the limit k -^ ±oo, we see that it is equal to 1, which 
gives (g]T|). Substituting this result back in (|]2B|), (|OB| ), we obtain ( ^ . 

Now to finish the proof it remains to obtain the local expression ( [4.9| ). To this end 
we first of all re- write the equation ( [4.15|) in terms of Sk and find the following difference 
equation: 

Sk+i - Sk = qk+irk+2 - qk-2rk-i- (4.27) 

As it stands, it does not allow a local solution. In order to achieve our goal, we have to 
once more switch on the equations of motion ( |4.8| ). Putting on the right-hand side of 

(P3) 

(gfc-i - qkSk+i){^ - qk+iTk+i) 
Qk+i — Qk+i — It 



(l + hD^'j [1 + hqkrk+2 + hqk-iVk+i - hqkVk+iSk+i 
[l + /i-Dfc ) /(l + hqk-iTk+i + hqk-2rk - hqk-iVkSj, 



Tk-i = rfc_i - h 
we obtain after some re- arrangements: 



(1 + hqkrk+2) 
(rfc+i - rkSk){l - qk-iTk-i) 



[Sk+i + hqk-.irk+2 



1 + hqkqk+irk+irk+2) 



[I + hqk-2rk) 



Sk + hqk-2rk+] 



Multiplying both sides by 



(1 + hqkrk+2) 
{qk+irk+2 - qk-2rk-i){^ + hqk-irk+i 

(1 -I- hqk-iqkrkrk+i) 



[1 + hqk-2qk-irk-irk) 
(1 + hqk-2rk) 



[1 + hqk-irk 



+ij 
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we arrive at the identity which may be read in the following way: the left-hand side of 
( [4. 91 ) does not depend on k. The value of this constant may be determined from the 
k — ^ ±00 limit. ■ 

Remark 1. It might be desirable to have the expressions for all entries of the matrix 
V^, in terms of (q, r) solely, and of the matrix V^ - in terms of (q, r). Notice that 
the quantities A),', Dl' already have the desired form. For the quantities -B^ , C^ it is 
possible to derive from ([4.2|) , (|4. 1|) the following nice expressions: 



p ^ , ~ , ~ 1 - hqkrk-2) ..^3. 

Pk - qk+iTk = (c + qkVk-i + qk~irk-2) jz Z ~ — ~ T' (4.28) 

(1 - hqkqk-irk-irk~2) 

D ~ ~ / , ~ I ~ ^ (-'■ ~ hqk+irk-i) , . 

Pk - qk-irk-2 = (c + qk+irk + qkrk-i) 7:^ r ^"^ r- (4-29) 

(1 - hqk+iqkrkrk-i) 

Analogously, for the quantities B^ , C^ it follows from ( [4.9|) , ( ^^) : 

c ^ , ~ , ~ (1 + hqk-2rk) . . „„s 

bk - qkrk+i = (c + qk-iTk + qk-2rk-i) ,^ , ,~ r, (4.30) 

(1 + hqk-2qk-irk-irk) 

c ~ ~ / , ~ , ~ ^ (1 + hqk-irk+i) . . 

Sk - qk-2rk-i = (c + qkTk+i + qk^iTk) ,-, , ,~ — r- (4.31) 

(1 + hqk-iqkr-krk+i) 

Unfortunately, analogous expressions for A^. , resp. D^ are complicated and non- 
elegant. 

Remark 2. Notice that from the Theorems 1,2 we can recover the maps from the 
previous paper [Q, together with the corresponding matrices Vk- Indeed, rescaling h \-^ 
h/c and then sending c —^ 00, we find from the Theorem 1 the map 

f {qk - qk)/h = -gfc+i(l - qkh), 
T,{-h) :{ ^ ^ (4.32) 

[ {rk - rk)/h = rk-i{l - qkVk) 

and the matrix 

/ l-h\^ + hqkh-i -hXqk 
V;«(q, ?,-/.)= |. (4.33) 

\ -h\rk-i 1 

and from the Theorem 2 the map 

{qk - qk)/h = -g/c-i(l - qkrk), 

T^ii-h) : { (4.34) 

{rk - rk)/h = rfc+i(l - qkTk) 
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along with the corresponding matrix 



.(-1) 



Vr'iSi,v.-h) 



1 h\ ^Qk-i 

hX'^Tk 1 + hX"'^ - hqk-iTk 



(4.35) 



Corollary 1. Consider the map 7^ ^(— /i;c) : (q, r) i-^ (q, r)- This map is defined by 
the equations of motion 



Tf\-h;c):l 



{Qk - qk)/h 
{h -rk)/h 



(gfc+2 - qk+iQk+i) (1 - qkrk) 

(1 + hqk+iTk-i) 
{rk-2 - Tk-iQk) (1 - qkTk) 



(4.36) 



where 



(1 + hqk+irk-i) 

(1 + hqkqk-irk-irk-2) (1 + hqk+iqkrkrk-i) 
(1 + hqkrk-2) (1 + hqk+iTk-i) 

-{qk+iTk + qkTk-i + qk-irk-2 + hqk+iqkqk-irkrk-irk-2) = c. 



(Qk + hq 



fc+l^fc-2 



It approximates the flow ( |2.18| ) and has a commutation representation ( |3.1j ) with the role 
of Vk played by the matrix 



W, 



(2) 



1 



/ 1 - hX^qkrk-i 



1 + hA 



(0) 



-hX^qk - hXB^f}^ \ 



(4.37) 



\^ -/iAVfc_i - hX&^> 1 + hA^^> + hX^A)f> - hX^ I 



1(2) E>(1) Wl) 4(0) 



(2) j.(l) ^(1) 



Here the quantities A\. , B^. , C\. , A\. are obtained from the quantities A\ , B)^ , C). , 
A\ by the change h to —h, q to q, and r to r. 

Corollary 2. Consider the map Tj^i—h^c) : (q, r) t-^ (q, ?)■ This map is defined by 
the equations of motion 



r-,\-h-c):{ 



where 



(Rk - hq 



'k~2^k+l 



(gfc-2 - qk-iRk) (1 - qkTk) 

(1 - hqk-ih+i) 

(^fc+2 - h+iRk+i) (1 - qkh) 

(1 - hqk-irk+i) 

1 - hqk-2qk-irk-ih) (1 - hqk-iqkhh+i] 
(1 - hqk-2rk) (1 - hqk-ih+i) 



{qk - qk)/h 
[rk - rk)/h 



(4.38) 
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■{qkrk+1 + Qk-iTk + qk-2rk-i - hqk-2qk-iqkrk-irkrk+i) = c. 



It approximates the flow ( |2.19| ) and has a commutation representation p.l| ) with the role 
of Vk played by the matrix 

I 1 - HDf - h\-^D[-^^ + h\'^ hX-^Qk-i + hX-^B'^^'-^ \ 

^ ^^k y hX-'h + hX-'ci-'^ l + hX-^k-ih J 

(4.39) 
Here the quantities D^. , B^. , Cj^ , Dk ^'"^ obtained from the quantities D)^ , 
Bj^ , Cj, , Df, by the change h to —h, q to q, and r to r. 

Proof. Since the both Corollaries 1,2 are proved analogously, we demonstrate again 
only the second one. The map TS2i—h;c) obviously allows the commutation represen- 
tation ( |3.1| ) with the matrix (V^ (q, r, —h)) playing the role of Vk. The Corollary 2 
will be proved, if we show that 

Wt'^ = il + hX-'){vt'\il,r,-h)y\ 

But from the expressions given in the Theorem 2 it is possible to derive that 

det Vt'\ci,r, h) = {l- hX'') (l + /^Z^f ; 

so that 

det Vt'\ii, r, -h) = (1 + hX-^) (l - hof^ 

which implies the above statement. ■ 

5 Relation to the Taha— Ablowitz schemes 

From the results of the previous Section it is easy to see that our maps are indeed some 
particular cases of the Taha-Ablowitz scheme. Namely: 

- The map T2{h] c) coincides with the Taha-Ablowitz difference scheme with the only 
nonzero parameters a^'^' = 1, a'-^-' = — c. 

- The map T^2{h]c) coincides with the Taha-Ablowitz difference scheme with the 
only nonzero parameters S^~^^ = —1, (5'-"^^ = c. 

- The map T^^^^—h; c) coincides with the Taha-Ablowitz difference scheme with the 
only nonzero parameters 5*^^-* = — 1, 5''^'* = c. 
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- The map T_2{—h; c) coincides with the Taha-Ablowitz difference scheme with the 
only nonzero parameters a^~^^ = 1, a*-"^-' = — c. 

These schemes, obviously, do not satisfy the condition (p.4|) (they certainly should not, 
because the underlying continuous time flows (|2.18|) , (|2.19|) do not allow the reduction 



]2[)). However, certain compositions of our maps, approximating the flow 

^2(t) o J'i(-ct) o J^_2(-t) o J^_i(ct), (5.1) 



do have the desired property ( |3.4|) , making them suitable for numerical integration of the 
space discretized mKdV equation ( |1.4|) . 

Namely, taking into account that the composition of maps corresponds to the mul- 
tiplication of the corresponding Vk matrices, we see immediately that the following two 
statements hold. 

Theorem 3. The map 

%{h-c)or^2{-h-c) (5.2) 

approximating the flow (|5.1| ), coincides with the Taha-Ablowitz difference scheme with the 
only nonzero parameters a^^' = 5^~^' = 1, a*-^^ = 5^^^' = —c, and allows the reduction 

Theorem 4. The map 

Tf\-h;c)oTS2\h;c) (5.3) 

approximating the flow (|5.1| ), coincides with the Taha-Ablowitz difference scheme with the 
only nonzero parameters a^"^^ = S^^'^ = —1, a*-~^^ = (5^^^ = c, and allows the reduction 

These two maps (taken by c = 2) are the first candidates to the role of valid difference 
schemes for numerical integration of the mKdV equation. Just as the scheme tested in 
0, they are integrable, i.e. they have infinitely many integrals of motion (in the rapidly 
decaying case; in the A^-periodic case the number of integrals reduces to N), and all the 
integrals are in involution with respect to the Posson bracket (p.8|) . Additionally, they have 
an advantage of locality. More precisely, each step of time integration by the corresponding 
numerical scheme consists of solving two local systems of nonlinear equations for the 
updates of q, r. We expect that the savings of the amount of computations due to 
locality will let these schemes to exceed the original Taha-Ablowitz scheme, even despite 
the fact that the latter is of the second order in h, while ( ^^) , (|5.3| ) are only of the first 
order in h. 

In this connection it should be mentioned that the maps ( |5.2|) , (|5.3| ) are suitable 
building blocks for applying the Ruth-Yoshida-Suzuki techniques []TU[, which in principle 
allow to construct the schemes of an arbitrary high order in h. The first step in this 
direction is given by the following statement. 
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Theorem 5. The composition 

T2{h/2; c) o T_2(-/i/2; c) o T-\-h/2; c) o TZ^{h/2- c) 



(5.4) 



approximates the flow ( ^.1|) up to the second order in h. It coincides with the Taha- 
Ahlowitz difference scheme with the parameters 



a 



a 



(2) 



^(-2) 



(4) ^ ^(-4) = _a(-4) = _ 
c{l-h/2) (_2) 

= 0, 



^(4) 



1 

2K' 



2A 

(0) 



<5(2) 



c(l + /i/2) 
2A 



a 



5(0) 



where 



A 



1 - /i'(l + c')/4 

In particular, this composition allows the reduction ( p.5| ). 

Proof. The map (|5.4]) , as a composition of two Taha-Ablowitz schemes, allows a 
commutation representation ( p.l|) with the matrix 14 whose A^; -D^ entries containing, 
in principle, terms with A-', j = 0, ±2, ±4, ±6, ±8, while the Bj., Ck entries contain the 
terms with j = ±1,±3,±5,±7. A careful inspection will convince that only the terms 
with |j| < 4 are actually present in this matrix. To this end we need the following two 
statements. 

Lemma 1. The entries of the matrix Vk for the map ( [5.2| ) have the following A- 
dependence: 



A{X) 
D{X) 
B{X) 
C(A) 



l + hX^ + hX^Af + hAf + hX-%^irk 
1 + hX^kh^i + /iDf) + /iA-^dI^') + hX-^ 
hX\k + hXB^^^ + hX-^^^-^^ + hX-%-1 



(1) 



,-ir^(-i) 



h\'7k-i + hXCl'' + hX-'C 



hX ^rk 



Lemma 2. 

dependence: 



The entries of the matrix Vk for the map ( ^l2|) have the following A- 



A{X) 
D{X) 
B{X) 
C{X) 



-2 A (-2) 
k 



Ak (l - hX%rk-i + /iA1°) + hX-^A 
Ak (l - hX" + hX^B^k^ + /iDf - hX-^ 



hX 



-4 



(1) 



,-1r(-i) 



Afc [hX% + hXB]:> + hX-'B 

Ak ihXhk^i + /lACi'^ + /lA-^Ci"'^ 



qk-iTk 

-3 



hX Qk^l 

+ hX-^h) 



15 



(The exact expressions for the unspecified functions play no role in the following 
reasonings; therefore we denoted carelessly different objects in the above statements by 
one and the same symbol, such as A^ ). The proof of these two lemmas are similar, so 
we give only the proof of the lemma 1. Let 

T_2{-h;c): (q,r)^(q,f), T2{h;c): (q,f)K^(q,r) 

It follows from the Theorems 1,2 that the the matrix Vk for the composition map ( p. 21 ) is 
equal to 

Vk = Vi'\^,r,h)vi-'\^,r,-h). 

The general A-dependence of this matrix follows directly from ([4 .31) , ([4.10|) , and we have 
only to prove the explicit expressions for B^ , C^ , A^ , D^ given in the formulation 
of the Lemma. The expressions C^ = ffc_i, C^ = r^ follow immediately, the rest is 
derived as follows: 

D?^ = Qkh-i {l ~ hDf^) + hh-iBi''^ = Qkh-i 
Bi''^ = qk-i {l - hAf^) + hBl}^ = q^-i 
aI"') = Qk-iTk {l - hA^^^) + hrkB^}^ = Qk-irk 

Here the first two expressions follow from ( [4.19|) (of course, with the change h to —h and 
q to q), and the last two expressions follow from an analogous "evolution equation for 
qk-i" for the map 7^. ■ 

Expressions given in the Lemmas 1,2 imply that in the product of the corresponding 
matrices Vk the following terms vanish identically: the terms with A^^, A^^ in the entries 
Ak, Dk, the terms with A''', A~^ in the entry B^, and the terms with A~^, A^ in the entry 
Cfc. After this the commutation representation ( p.l| ) shows that the last obstacles, namely 
the term with A^ in Bk and the term with A~^ in Ck, also must vanish. So, the matrix Vk 



corresponding to the composition ( p.4| ) has the desired A-dependence. Now it remains to 



observe that the limit values of its entries by A; — > ±oo are: 

lim Ak{X) = (l- hX-^/2 + h\-^c/2) (l + h\^ /2 - h\^c/2 



= A (l + /lA^a(^) + /iA2a(2) + /,a-2«(-2) + /.a-^^C-^) 
lim Dk{X) = (l - /iA72 + hX^c/2) (l + hX-^/2 - hX~^c/2 

fe— >±oo ^ ^ ^ 

= A (l + /lA^^W + h\H^^^ + h\-H^-^^ + /iA-^(5(-^)) 
with the values a^^\ 6^^^ given in the Theorem. ■ 
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6 Conclusion 

In the present paper we re-considered the discretizations of the modified Korteweg-de 
Vries equation due to Taha and Ablowitz. We demonstrated that by some choice of 
parameters their highly non-local scheme may be factorized into the product of much 
more simple (local) ones, each of them approximating a more simple and fundamental 
fiow of the Ablowitz-Ladik hierarchy. These local schemes may be studied exhaustively. 
In particular, by the same change of variables as in the previous paper [Q, we can establish 
a relation to the pair of "second" discrete time fiows of the relativistic Toda hierarchy, 
which gives a simple way to determine interpolating Hamiltonian fiows for our maps 
and to solve them in terms of a factorization problem in a loop group (the relation 
between the Ablowitz-Ladik and the relativistic Toda hierarchies is due to |]TT|]). We 
guess that also in the practical computations our variant of the difference scheme will 
exceed considerably the old one. It would be interesting and important to carry out the 
corresponding numerical experiments. 

The research of the author is financially supported by the DFG (Deutsche Forschungs- 
gemeinshaft). 
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7 Appendix: 

difference scheme by Taha and Ablowitz 



owing from (|3l^)-(|3l9|), read: 


.(4) 


-4^' 


= 







4 (-4) 


■ 4-^' 


= 


*c<-> . 


- nBl;^ 


4(2) 
^k+1 


-4^' 


= 


*c«' - 


r.BS, 


4{-2) 


■ 4-^' 


= 


*c<-'> . 


- nBl;^ 


^k+i 


-4°' 


= 


*cl" - 


^k^k+1 


n(-4) _ 
^k+i 


-4-" 


^ 







r)(4) 


-i^i" 


— 


?.bF' - 


n r(3) 


n(-2) _ 
^k+i 


-D[r-'' 


— 


hBt^' 


n r(-3) 


r)(2) 


'D^' 


— 


hBl^> - 


IkC^^l 


flSi 


-Df 


— 


r.Br 


- qkCk+i 


Bf 


= qtA 




- *Bf 





[11] S.Kharchev, A.Mironov, A.Zhedanov. Faces of relativistic Toda chain. - |hep- 
I th/9606T4^ . 



(7.1) 

(7.2) 
(7.3) 
(7.4) 
(7.5) 

(7.6) 
(7.7) 
(7.8) 
(7.9) 
(7.10) 

(7.11) 
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B 



(-3) 
fc+1 






1k^k+l 



B 



(1) 



B 



(-1) 
fc+i 



gyt^i+i - qkD\ 

Qk^k 



^(3) 



„ 4(-2) , r(-3) 






(-3) 






(3) 
fc+1 



Qk -Qk 

Tk-Tk 



(1) 

fc+1 



'fc-'^fc+i 

r, 4(4) _ r)(4) 

rfc4'^ - r.Dfl 



C, 



(-3) 
fc+1 



c, 



(3) 



h(n 4(°) ;? nW -u r(i) r(-i 
/i(^gfc/ifc+i - ik^k + ^fc+i - ^fc 



)(0) 



(0) 



.(-1) 



HTkD^kU-rkAT + Cr^'^-C, 



^(1) 



7.12) 
7.13) 
7.14) 

7.15) 
7.16) 
7.17) 
7.18) 

7.19) 
7.20) 



Taha and Ablowitz interpret the first 18 equations ( |7.1| )-( jrTE| ) as the defining relations 
for the 18 coefficients A^ -D^\ and the last two equations ( |7.19| ), ( |7.20D as the evolution 
equations for g^, fk- This leads to a highly nonlocal scheme. Let us indicate how do Taha 
and Ablowitz solve the first 18 equations (with some amendments). 

First of all, (^, Q imply 



A 



(4) 



a 



(4) 



D 



(-4) 



^(-4) 



(7.21) 



where a'^^\ 5^ ^^ are some constants (and the same holds for other a'^^\ 5^^^ appearing 
further on). Substituting ( fTATD , ( |7A6| ) into (^), and (fTT^ ), {^A§ into (|73), we get 
upon use of ( [7.21| ) the following difference equations: 



l(4) 



D):_1,{1 - q,r,) 



Dl^\l-q,h) 



l(-4) 



4+Ai-gfcr,) 



Ai-'\l-q,h) 



Their solution may be put into the form 



D 



(4) 



5(^)A, 



A 



(-4) 



a 



(-^^Afc, 



(7.22) 



(7.23) 



where 



fc-i 



A, 



j=-oo 



(Ij^j 



Wi 



(7.24) 



is for all k of order 1 + 0{h) and tends to 1 as A; —i> ±00 (due to the infinite- dimensional 
analog of the integral of motion 
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The back substitution (g), Ull^ into ([77[T| ), ([71^) , ([71^ ), (^) results in 



5f ) = a(^)g, - 6'^'^qkAk, (7.25) 

^M) ^ ^(-4)^^ _ «(-4)^^A,, (7.26) 

5(-3) ^ 5(-4)~^_^ _ tt(-4)g,_^A,, (7.27) 

Cf = «(^)Ffc_i - 5(^Vfc_iA,. (7.28) 



Substituting the last formulas into ( [OD , ( |f^ , we arrive at the following difference equa- 
tions: 

^1+1 - ^k^ = a'-^Kqkrk-i - qk+iTk) + ^^'^Kqk+iTk^k+i - Qkrk-i^k), (7.29) 

D^+i - Dir^^ = S^~^\qk-irk - qkTk+i) + a^-'^\qkrk+iK+i - qk-ihAk)- (7.30) 

Their solutions are given by 

Af = «(2) _ a^^\qkrk-i + $fc) + ^^^^g^r.^iA^, (7.31) 

^(-2) = ^(-2) _ ^(-4) (^^^^^^ + v]/,) + a(-4)g,_,FfcAfc, (7.32) 

where the quantites 

fc-i fc-i 

'^fc= II (^i'^i-i-^i^j-i)) ^fc= II (?i-i^i - ^^i-i^i) (7-33) 

j = — 00 j = — OD 

are of order 0{h) and vanish by A; ^ ±oo (the statement for k —>■ +oo follows from the 
fact that the infinite dimensional analogs of the quantities if±i(q, r) from (|2.11|) are the 
integrals of motion) . 

To proceed further, we substitute ( |7A3|) , ( [TlSD into (|7]|). Using (^^25]), (3), and 



( [7.29| ), we derive from (|7.9|) the following difference equation: 

{D^l^ - a('')gfc+i^fc)(l - qkTk) - {Dk^ - a^^^gfcrfc_i)(l - qkh) = 

-S'^^^ (qk+ihil - qkrk)Ak+i - qkrk-i{l - qkrk)Ak)- 
Analogously, substituting ( fTTD , (ITTtD into Q and using ( [^271) , ( f7^ , ( [^301) , we 



obtain: 

(4;? - (5(-^)gfcr,+i)(l - qkn) - (4"'^ - S^-'^qk^irk){l - qkh) 
if("'^)(gfcffc+i(l - qkrk)Ak+i - qk-irk{l - qkrk)A]}j. 
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-a' 



A little trick (apparently missed in 0), based on the formula {l—qkrk)Ak^i = {l—q^rkiKk, 
allows to solve these difference equations as 

Df = 5(2) A, + «(4)g,ffc_i - 5^'^\qurk-i - $fc)Afc, (7.34) 

4"') = a(-2) A, + 5^-%k-irk - a^'^^ {qk^.h - ^fc) A^. (7.35) 

Now we have determined all the quantities from the right-hand sides of ( |7.13| ) , ( [7. 141 ) , 
( [r.lSp , ( [7.171) . Collecting all the results, we find: 



4'^ = «(2)g,-<5(2)g,A, 

+a^^^(gfc+i - Qkiqk+irk + QkTk-i + ^k)] 
-d^'^^ijk+i - qkiqk+ih + qkh^i - ^k))^k (7.36) 

Bt^ = 6^~'^qk-i-a^''kk^iAk 

+S^~'^^{qk-2 - qk-iiqk-iTk + qk-2rk-i + ^fc)) 

-a^^^^ {qk^2 - qk-i{qk-2rk~i + qk-ih - ^fe)) A^ (7.37) 

Cj^^ = «(^)r,_i - <5(2)r,_iA, 

+a^'^^{rk-2 - Tk-iiqkTk-i + qk-irk-2 + ^k)) 

-6^'^Urk-2 - rk-i{qk~irk-2 + qkh-i - $fc))Afc (7.38) 



Ct'^ = 5^-\k-a^-'^hAk 

+(5^""^) (rfc+i - VkiqkTk+i + qk^iTk + ^k)\ 

-a^"^) (ffc+i - hiqkh+i + qk-iTk - ^fc)) Afc (7.39) 



Finally, these expressions, being substituted in (|7.5|) , ( |7.1CI|) , result in difference equa- 
tions whose solutions may be represented as 

4°) = a(0) - a'^'\qkrk^i + <^k) + S^''>qkrk-,Ak 

-a^'^\qk+irk-i + qkrk-2 - qk+iqkhrk-i - qkqkh-irk^i - qkqk^irk~irk-2 + Sfc) 
+ 6^^^^ (qk+iVk-i + qkrk-2 - qkrk-iiqk+ih + qkTk-i + qk-irk-2 - ^k))Ak (7.40) 

Df = 5(°) - 5(-2)(gfc„ir, + vj/fc) + a^-'^k-ihAk 
-5^~'^\qk-irk+i + qk-2rk - qkqk-iTk+iTk - qk-iqk-ihrk - qk-iqk-2rkrk-i + ^k) 
+ a^""^) (gfc_iffc+i + qk-2rk - qk-ihiqkh+i + qk-ih + qk-2rk-i - ^fc)) A^ (7.41) 
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Here 

fc-i 

j=-oo 

fc-1 

j=-co 

These quantities, just as $fc, \l/fc, are of order 0{h) and vanish in both hmits k — > ±cx). 
The see that this is true for k — > +cxd, one has to use the partial summation. For example, 
for Sfc one arrives at the following expression: 

Sfc = qkrk-i^k+2 - -{qkrk-i<Pk-i + qkrk-i<Pk-i - 2qkrk-i4>k-i) 

j=-oo 

fc-1 1 

- E {wj-2 - i^wj-i^qj+i^j + Wj-i + qj-1^3-2)). 



where 



j=-oo 

fc-1 



j=-oo 

The obtained expression for S^ clearly vanishes at A; ^ +00 due to the fact that an 
infinite-dimensional analog of (p.l2|) is an integral of motion. 

Now the Taha-Ablowitz difference scheme for the mKdV is given by ( |7.19| ) , ( |7.20| ) with 



the expressions (|7.36|) - (|7.41|) . Due to the appearance of the quantities A^, <l>fc, ^fc, S^, Vtk 
this scheme is not only implicit, but also highly non-local. More precisely, they consider 
only the reduced case r = ±q, assuming 5'^^^ = a'^~^\ In this reduction A)!' = DjT , 
Bf^ = ±Cjr , and, of course, $fc = \E'fc, S^ = fi^, so that the number of the nonlocal 
ingredients is reduced by two. 
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